Import model objects
load('/Users/hp2500/Google Drive/STUDY/Columbia/Research/Corona/Data/GER/ger_list_results_fixed_window.RData')
load('/Users/hp2500/Google Drive/STUDY/Columbia/Research/Corona/Data/US/us_list_results_fixed_window.RData')
Define Functions
list_iterater <- function(models, test) {
for(i in models){
for(j in i){
if(test == 'qq'){j %>% plot(2)}
if(test == 'ks'){j %>% resid() %>% ks.test(y=pnorm) %>% print()}
if(test == 'bp'){j %>% bptest() %>% print()}
if(test == 'ph'){j %>% cox.zph() %>% print()}
}
}
}
Assumptions GER COVID-19 onsets (proportional hazards)
list_iterater(ger_list_results$ger_cox_prev_onset, test = 'ph')
chisq df p
pers 6.21 1 0.013
GLOBAL 6.21 1 0.013
chisq df p
pers 4.007 1 0.045
age 21.844 1 3.0e-06
male 0.454 1 0.500
conservative 15.769 1 7.2e-05
GLOBAL 25.231 4 4.5e-05
chisq df p
pers 1.30051 1 0.25
academics 0.10832 1 0.74
medinc 0.00668 1 0.93
manufact 0.07463 1 0.78
GLOBAL 1.38751 4 0.85
chisq df p
pers 6.89 1 0.0087
airport_dist 7.94 1 0.0048
tourism 1.65 1 0.1994
healthcare 6.31 1 0.0120
popdens 6.31 1 0.0120
GLOBAL 17.20 5 0.0041
chisq df p
pers 1.05049 1 0.30539
age 11.81505 1 0.00059
male 4.76514 1 0.02904
conservative 4.75085 1 0.02928
academics 0.11884 1 0.73030
medinc 2.27744 1 0.13127
manufact 0.49479 1 0.48180
airport_dist 0.69603 1 0.40412
tourism 0.07452 1 0.78486
healthcare 0.00138 1 0.97032
popdens 1.36622 1 0.24246
GLOBAL 17.21722 11 0.10161
chisq df p
pers 2.51 1 0.11
GLOBAL 2.51 1 0.11
chisq df p
pers 7.888 1 0.00498
age 19.642 1 9.3e-06
male 0.658 1 0.41721
conservative 14.163 1 0.00017
GLOBAL 23.934 4 8.2e-05
chisq df p
pers 1.6479 1 0.20
academics 0.0923 1 0.76
medinc 0.0737 1 0.79
manufact 0.0601 1 0.81
GLOBAL 1.9272 4 0.75
chisq df p
pers 3.23 1 0.0723
airport_dist 8.82 1 0.0030
tourism 1.55 1 0.2126
healthcare 6.01 1 0.0142
popdens 6.41 1 0.0114
GLOBAL 16.26 5 0.0061
chisq df p
pers 3.45e+00 1 0.06333
age 1.28e+01 1 0.00034
male 4.89e+00 1 0.02701
conservative 5.00e+00 1 0.02530
academics 8.68e-02 1 0.76827
medinc 2.30e+00 1 0.12979
manufact 6.41e-01 1 0.42342
airport_dist 8.22e-01 1 0.36449
tourism 6.81e-02 1 0.79418
healthcare 9.25e-04 1 0.97574
popdens 1.38e+00 1 0.23943
GLOBAL 1.87e+01 11 0.06668
chisq df p
pers 10.2 1 0.0014
GLOBAL 10.2 1 0.0014
chisq df p
pers 9.075 1 0.0026
age 20.874 1 4.9e-06
male 0.256 1 0.6128
conservative 15.413 1 8.6e-05
GLOBAL 29.120 4 7.4e-06
chisq df p
pers 3.3334 1 0.068
academics 0.3435 1 0.558
medinc 0.0519 1 0.820
manufact 0.0950 1 0.758
GLOBAL 3.3556 4 0.500
chisq df p
pers 8.30 1 0.0040
airport_dist 7.50 1 0.0062
tourism 1.70 1 0.1924
healthcare 7.55 1 0.0060
popdens 6.36 1 0.0117
GLOBAL 21.03 5 0.0008
chisq df p
pers 4.09e+00 1 0.04323
age 1.17e+01 1 0.00064
male 4.38e+00 1 0.03633
conservative 4.86e+00 1 0.02754
academics 2.09e-01 1 0.64747
medinc 2.49e+00 1 0.11445
manufact 4.16e-01 1 0.51913
airport_dist 7.69e-01 1 0.38049
tourism 8.12e-02 1 0.77565
healthcare 5.23e-04 1 0.98176
popdens 1.31e+00 1 0.25286
GLOBAL 1.83e+01 11 0.07556
chisq df p
pers 7.14 1 0.0076
GLOBAL 7.14 1 0.0076
chisq df p
pers 3.46 1 0.0628
age 21.13 1 4.3e-06
male 0.44 1 0.5072
conservative 15.05 1 0.0001
GLOBAL 28.94 4 8.1e-06
chisq df p
pers 2.9859 1 0.084
academics 0.0492 1 0.825
medinc 0.0130 1 0.909
manufact 0.0662 1 0.797
GLOBAL 3.1061 4 0.540
chisq df p
pers 7.98 1 0.00474
airport_dist 8.47 1 0.00362
tourism 1.55 1 0.21352
healthcare 5.87 1 0.01537
popdens 6.24 1 0.01247
GLOBAL 20.74 5 0.00091
chisq df p
pers 2.34952 1 0.1253
age 12.10476 1 0.0005
male 5.08073 1 0.0242
conservative 4.84600 1 0.0277
academics 0.03360 1 0.8546
medinc 2.16525 1 0.1412
manufact 0.63890 1 0.4241
airport_dist 0.64944 1 0.4203
tourism 0.06193 1 0.8035
healthcare 0.00924 1 0.9234
popdens 1.17671 1 0.2780
GLOBAL 21.93744 11 0.0249
chisq df p
pers 5.91 1 0.015
GLOBAL 5.91 1 0.015
chisq df p
pers 4.570 1 0.033
age 19.880 1 8.2e-06
male 0.332 1 0.564
conservative 15.516 1 8.2e-05
GLOBAL 26.962 4 2.0e-05
chisq df p
pers 0.507 1 0.48
academics 0.202 1 0.65
medinc 0.128 1 0.72
manufact 0.280 1 0.60
GLOBAL 1.003 4 0.91
chisq df p
pers 5.90 1 0.01512
airport_dist 7.85 1 0.00509
tourism 1.40 1 0.23743
healthcare 5.39 1 0.02023
popdens 7.23 1 0.00716
GLOBAL 20.61 5 0.00096
chisq df p
pers 0.7484 1 0.38699
age 10.8876 1 0.00097
male 4.0270 1 0.04478
conservative 5.2602 1 0.02182
academics 0.2251 1 0.63517
medinc 1.1519 1 0.28314
manufact 0.1375 1 0.71080
airport_dist 0.8884 1 0.34590
tourism 0.0349 1 0.85181
healthcare 0.0493 1 0.82426
popdens 1.4904 1 0.22216
GLOBAL 15.5926 11 0.15694
Assumptions US COVID-19 onsets (proportional hazards)
list_iterater(us_list_results$us_cox_prev_onset, test = 'ph')
chisq df p
pers 133 1 <2e-16
GLOBAL 133 1 <2e-16
chisq df p
pers 106.580 1 < 2e-16
age 0.453 1 0.5
male 21.391 1 3.7e-06
conservative 102.647 1 < 2e-16
GLOBAL 172.622 4 < 2e-16
chisq df p
pers 96.8 1 < 2e-16
academics 144.2 1 < 2e-16
medinc 47.3 1 5.9e-12
manufact 61.3 1 5.0e-15
GLOBAL 182.7 4 < 2e-16
chisq df p
pers 92.07 1 < 2e-16
airport_dist 8.83 1 0.003
tourism 15.46 1 8.4e-05
healthcare 37.50 1 9.1e-10
popdens 3.87 1 0.049
GLOBAL 122.81 5 < 2e-16
chisq df p
pers 71.8062 1 < 2e-16
age 0.0164 1 0.89806
male 21.4651 1 3.6e-06
conservative 70.0537 1 < 2e-16
academics 110.6622 1 < 2e-16
medinc 32.3876 1 1.3e-08
manufact 46.1464 1 1.1e-11
airport_dist 2.0257 1 0.15466
tourism 13.2975 1 0.00027
healthcare 25.6451 1 4.1e-07
popdens 32.9347 1 9.5e-09
GLOBAL 173.6175 11 < 2e-16
chisq df p
pers 0.729 1 0.39
GLOBAL 0.729 1 0.39
chisq df p
pers 2.56 1 0.11
age 1.30 1 0.25
male 20.65 1 5.5e-06
conservative 107.90 1 < 2e-16
GLOBAL 124.50 4 < 2e-16
chisq df p
pers 0.112 1 0.74
academics 153.781 1 < 2e-16
medinc 53.480 1 2.6e-13
manufact 64.290 1 1.1e-15
GLOBAL 170.151 4 < 2e-16
chisq df p
pers 0.425 1 0.5143
airport_dist 7.913 1 0.0049
tourism 15.657 1 7.6e-05
healthcare 41.027 1 1.5e-10
popdens 1.702 1 0.1921
GLOBAL 59.516 5 1.5e-11
chisq df p
pers 2.14e-01 1 0.64358
age 4.78e-03 1 0.94491
male 1.97e+01 1 9.2e-06
conservative 7.15e+01 1 < 2e-16
academics 1.17e+02 1 < 2e-16
medinc 3.90e+01 1 4.2e-10
manufact 4.88e+01 1 2.8e-12
airport_dist 1.27e+00 1 0.26002
tourism 1.39e+01 1 0.00019
healthcare 2.87e+01 1 8.3e-08
popdens 3.04e+01 1 3.5e-08
GLOBAL 1.71e+02 11 < 2e-16
chisq df p
pers 1.91 1 0.17
GLOBAL 1.91 1 0.17
chisq df p
pers 0.913 1 0.34
age 1.251 1 0.26
male 22.838 1 1.8e-06
conservative 111.394 1 < 2e-16
GLOBAL 124.947 4 < 2e-16
chisq df p
pers 0.964 1 0.33
academics 148.095 1 < 2e-16
medinc 51.959 1 5.7e-13
manufact 62.359 1 2.9e-15
GLOBAL 162.064 4 < 2e-16
chisq df p
pers 0.831 1 0.3619
airport_dist 8.887 1 0.0029
tourism 13.795 1 0.0002
healthcare 39.805 1 2.8e-10
popdens 3.149 1 0.0760
GLOBAL 58.645 5 2.3e-11
chisq df p
pers 2.52e-01 1 0.61546
age 4.88e-03 1 0.94429
male 2.27e+01 1 1.9e-06
conservative 7.45e+01 1 < 2e-16
academics 1.14e+02 1 < 2e-16
medinc 3.65e+01 1 1.6e-09
manufact 4.77e+01 1 4.9e-12
airport_dist 1.94e+00 1 0.16421
tourism 1.26e+01 1 0.00038
healthcare 2.73e+01 1 1.7e-07
popdens 3.49e+01 1 3.6e-09
GLOBAL 1.69e+02 11 < 2e-16
chisq df p
pers 0.402 1 0.53
GLOBAL 0.402 1 0.53
chisq df p
pers 1.03 1 0.31
age 1.37 1 0.24
male 21.87 1 2.9e-06
conservative 111.50 1 < 2e-16
GLOBAL 134.61 4 < 2e-16
chisq df p
pers 0.603 1 0.44
academics 153.090 1 < 2e-16
medinc 54.952 1 1.2e-13
manufact 68.508 1 < 2e-16
GLOBAL 172.867 4 < 2e-16
chisq df p
pers 0.0332 1 0.8554
airport_dist 8.7075 1 0.0032
tourism 15.3470 1 8.9e-05
healthcare 42.5992 1 6.7e-11
popdens 3.5426 1 0.0598
GLOBAL 61.8373 5 5.1e-12
chisq df p
pers 3.81e-03 1 0.95076
age 1.66e-04 1 0.98973
male 2.11e+01 1 4.4e-06
conservative 7.67e+01 1 < 2e-16
academics 1.19e+02 1 < 2e-16
medinc 3.83e+01 1 6.1e-10
manufact 5.16e+01 1 6.7e-13
airport_dist 1.77e+00 1 0.18394
tourism 1.36e+01 1 0.00023
healthcare 3.01e+01 1 4.2e-08
popdens 3.77e+01 1 8.1e-10
GLOBAL 1.76e+02 11 < 2e-16
chisq df p
pers 49.7 1 1.8e-12
GLOBAL 49.7 1 1.8e-12
chisq df p
pers 35.610 1 2.4e-09
age 0.744 1 0.39
male 21.775 1 3.1e-06
conservative 109.483 1 < 2e-16
GLOBAL 133.700 4 < 2e-16
chisq df p
pers 47.0 1 7.2e-12
academics 158.3 1 < 2e-16
medinc 58.3 1 2.2e-14
manufact 62.3 1 2.9e-15
GLOBAL 177.2 4 < 2e-16
chisq df p
pers 33.48 1 7.2e-09
airport_dist 6.46 1 0.01102
tourism 12.01 1 0.00053
healthcare 41.31 1 1.3e-10
popdens 3.34 1 0.06744
GLOBAL 81.47 5 4.1e-16
chisq df p
pers 25.1056 1 5.4e-07
age 0.0324 1 0.85712
male 21.4504 1 3.6e-06
conservative 71.5756 1 < 2e-16
academics 119.2835 1 < 2e-16
medinc 40.8184 1 1.7e-10
manufact 47.3892 1 5.8e-12
airport_dist 1.6077 1 0.20481
tourism 11.5992 1 0.00066
healthcare 28.6560 1 8.6e-08
popdens 30.7441 1 2.9e-08
GLOBAL 176.2965 11 < 2e-16
Assumptions GER COVID-19 growth rates (normality of residuals)
list_iterater(ger_list_results$ger_lm_prev_slope, test = 'qq')

























list_iterater(ger_list_results$ger_lm_prev_slope, test = 'bp')
studentized Breusch-Pagan test
data: .
BP = 0.20014, df = 1, p-value = 0.6546
studentized Breusch-Pagan test
data: .
BP = 1.7565, df = 4, p-value = 0.7804
studentized Breusch-Pagan test
data: .
BP = 8.8926, df = 4, p-value = 0.06384
studentized Breusch-Pagan test
data: .
BP = 12.765, df = 5, p-value = 0.02568
studentized Breusch-Pagan test
data: .
BP = 19.536, df = 11, p-value = 0.05213
studentized Breusch-Pagan test
data: .
BP = 0.014025, df = 1, p-value = 0.9057
studentized Breusch-Pagan test
data: .
BP = 2.0882, df = 4, p-value = 0.7195
studentized Breusch-Pagan test
data: .
BP = 8.2989, df = 4, p-value = 0.08122
studentized Breusch-Pagan test
data: .
BP = 10.635, df = 5, p-value = 0.05912
studentized Breusch-Pagan test
data: .
BP = 18.827, df = 11, p-value = 0.06426
studentized Breusch-Pagan test
data: .
BP = 0.9348, df = 1, p-value = 0.3336
studentized Breusch-Pagan test
data: .
BP = 2.8075, df = 4, p-value = 0.5905
studentized Breusch-Pagan test
data: .
BP = 8.2707, df = 4, p-value = 0.08215
studentized Breusch-Pagan test
data: .
BP = 10.975, df = 5, p-value = 0.05188
studentized Breusch-Pagan test
data: .
BP = 18.978, df = 11, p-value = 0.06149
studentized Breusch-Pagan test
data: .
BP = 1.9621, df = 1, p-value = 0.1613
studentized Breusch-Pagan test
data: .
BP = 2.3503, df = 4, p-value = 0.6716
studentized Breusch-Pagan test
data: .
BP = 8.6682, df = 4, p-value = 0.06995
studentized Breusch-Pagan test
data: .
BP = 17.72, df = 5, p-value = 0.003318
studentized Breusch-Pagan test
data: .
BP = 19.885, df = 11, p-value = 0.04694
studentized Breusch-Pagan test
data: .
BP = 0.77295, df = 1, p-value = 0.3793
studentized Breusch-Pagan test
data: .
BP = 2.6459, df = 4, p-value = 0.6187
studentized Breusch-Pagan test
data: .
BP = 9.1115, df = 4, p-value = 0.05837
studentized Breusch-Pagan test
data: .
BP = 11.656, df = 5, p-value = 0.03981
studentized Breusch-Pagan test
data: .
BP = 19.53, df = 11, p-value = 0.05222
Assumptions US COVID-19 growth rates (normality of residuals)
list_iterater(us_list_results$us_lm_prev_slope, test = 'qq')

























list_iterater(us_list_results$us_lm_prev_slope, test = 'bp')
studentized Breusch-Pagan test
data: .
BP = 25.881, df = 1, p-value = 3.63e-07
studentized Breusch-Pagan test
data: .
BP = 58.996, df = 4, p-value = 4.715e-12
studentized Breusch-Pagan test
data: .
BP = 27.503, df = 4, p-value = 1.573e-05
studentized Breusch-Pagan test
data: .
BP = 88.757, df = 5, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 121.41, df = 11, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 0.45532, df = 1, p-value = 0.4998
studentized Breusch-Pagan test
data: .
BP = 58.951, df = 4, p-value = 4.819e-12
studentized Breusch-Pagan test
data: .
BP = 35.857, df = 4, p-value = 3.096e-07
studentized Breusch-Pagan test
data: .
BP = 93.365, df = 5, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 127.65, df = 11, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 3.7171, df = 1, p-value = 0.05386
studentized Breusch-Pagan test
data: .
BP = 59.824, df = 4, p-value = 3.158e-12
studentized Breusch-Pagan test
data: .
BP = 21.132, df = 4, p-value = 0.0002982
studentized Breusch-Pagan test
data: .
BP = 84.501, df = 5, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 119.57, df = 11, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 7.0583, df = 1, p-value = 0.00789
studentized Breusch-Pagan test
data: .
BP = 65.59, df = 4, p-value = 1.933e-13
studentized Breusch-Pagan test
data: .
BP = 57.311, df = 4, p-value = 1.064e-11
studentized Breusch-Pagan test
data: .
BP = 108.23, df = 5, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 147.26, df = 11, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 12.544, df = 1, p-value = 0.0003975
studentized Breusch-Pagan test
data: .
BP = 58.66, df = 4, p-value = 5.545e-12
studentized Breusch-Pagan test
data: .
BP = 31.917, df = 4, p-value = 1.989e-06
studentized Breusch-Pagan test
data: .
BP = 94.189, df = 5, p-value < 2.2e-16
studentized Breusch-Pagan test
data: .
BP = 123.74, df = 11, p-value < 2.2e-16
Assumptions GER socdist onsets
list_iterater(ger_list_results$ger_cox_socdist_cpt, test = 'ph')
chisq df p
pers 0.187 1 0.67
GLOBAL 0.187 1 0.67
chisq df p
pers 0.589 1 0.443
age 1.510 1 0.219
male 5.639 1 0.018
conservative 1.333 1 0.248
GLOBAL 8.287 4 0.082
chisq df p
pers 0.07128 1 0.789
academics 0.00457 1 0.946
medinc 6.30499 1 0.012
manufact 5.40654 1 0.020
GLOBAL 9.07479 4 0.059
chisq df p
pers 0.0093 1 0.923
airport_dist 3.7655 1 0.052
tourism 0.4039 1 0.525
healthcare 0.9540 1 0.329
popdens 0.0948 1 0.758
GLOBAL 5.0000 5 0.416
chisq df p
pers 0.0439 1 0.834
age 3.8044 1 0.051
male 6.2038 1 0.013
conservative 2.3823 1 0.123
academics 0.0481 1 0.826
medinc 3.3040 1 0.069
manufact 2.3737 1 0.123
airport_dist 6.1123 1 0.013
tourism 1.5087 1 0.219
healthcare 1.5767 1 0.209
popdens 0.0263 1 0.871
onset_prev 2.8203 1 0.093
slope_prev 4.1380 1 0.042
GLOBAL 17.3136 13 0.185
chisq df p
pers 0.284 1 0.59
GLOBAL 0.284 1 0.59
chisq df p
pers 1.15 1 0.284
age 1.34 1 0.247
male 5.15 1 0.023
conservative 1.12 1 0.290
GLOBAL 7.65 4 0.105
chisq df p
pers 0.5249 1 0.469
academics 0.0037 1 0.952
medinc 5.4471 1 0.020
manufact 5.4390 1 0.020
GLOBAL 7.8952 4 0.095
chisq df p
pers 0.1079 1 0.743
airport_dist 4.6196 1 0.032
tourism 0.4192 1 0.517
healthcare 1.2591 1 0.262
popdens 0.0509 1 0.822
GLOBAL 6.7462 5 0.240
chisq df p
pers 0.5092 1 0.4755
age 3.8442 1 0.0499
male 5.7297 1 0.0167
conservative 2.3451 1 0.1257
academics 0.0958 1 0.7570
medinc 3.2215 1 0.0727
manufact 2.2607 1 0.1327
airport_dist 6.8590 1 0.0088
tourism 1.4968 1 0.2212
healthcare 1.6372 1 0.2007
popdens 0.0114 1 0.9150
onset_prev 3.0881 1 0.0789
slope_prev 4.5326 1 0.0333
GLOBAL 17.6092 13 0.1729
chisq df p
pers 1.04 1 0.31
GLOBAL 1.04 1 0.31
chisq df p
pers 1.52 1 0.218
age 1.18 1 0.276
male 4.82 1 0.028
conservative 1.05 1 0.307
GLOBAL 7.57 4 0.109
chisq df p
pers 0.9332 1 0.334
academics 0.0254 1 0.873
medinc 5.5409 1 0.019
manufact 4.7740 1 0.029
GLOBAL 8.6942 4 0.069
chisq df p
pers 1.232 1 0.27
airport_dist 4.230 1 0.04
tourism 0.475 1 0.49
healthcare 1.037 1 0.31
popdens 0.103 1 0.75
GLOBAL 5.784 5 0.33
chisq df p
pers 0.8818 1 0.3477
age 3.7021 1 0.0543
male 5.4792 1 0.0192
conservative 2.2946 1 0.1298
academics 0.1120 1 0.7379
medinc 2.5620 1 0.1095
manufact 1.7891 1 0.1810
airport_dist 7.1031 1 0.0077
tourism 1.7503 1 0.1858
healthcare 1.6702 1 0.1962
popdens 0.0123 1 0.9115
onset_prev 2.9767 1 0.0845
slope_prev 4.3407 1 0.0372
GLOBAL 17.9668 13 0.1588
chisq df p
pers 0.667 1 0.41
GLOBAL 0.667 1 0.41
chisq df p
pers 1.50 1 0.221
age 1.30 1 0.255
male 4.96 1 0.026
conservative 1.09 1 0.295
GLOBAL 8.17 4 0.086
chisq df p
pers 0.96909 1 0.325
academics 0.00695 1 0.934
medinc 5.51116 1 0.019
manufact 5.08746 1 0.024
GLOBAL 7.82380 4 0.098
chisq df p
pers 0.4477 1 0.503
airport_dist 4.4029 1 0.036
tourism 0.4347 1 0.510
healthcare 1.0666 1 0.302
popdens 0.0851 1 0.770
GLOBAL 5.8358 5 0.323
chisq df p
pers 1.3945 1 0.2377
age 3.7995 1 0.0513
male 5.6669 1 0.0173
conservative 2.3341 1 0.1266
academics 0.0935 1 0.7598
medinc 3.2418 1 0.0718
manufact 2.2310 1 0.1353
airport_dist 6.8871 1 0.0087
tourism 1.5300 1 0.2161
healthcare 1.6230 1 0.2027
popdens 0.0108 1 0.9173
onset_prev 2.9882 1 0.0839
slope_prev 4.5077 1 0.0337
GLOBAL 17.5141 13 0.1769
chisq df p
pers 0.66 1 0.42
GLOBAL 0.66 1 0.42
chisq df p
pers 0.268 1 0.605
age 1.388 1 0.239
male 4.860 1 0.027
conservative 1.229 1 0.268
GLOBAL 6.221 4 0.183
chisq df p
pers 0.91068 1 0.340
academics 0.00965 1 0.922
medinc 5.97668 1 0.014
manufact 5.24467 1 0.022
GLOBAL 7.79339 4 0.099
chisq df p
pers 0.6571 1 0.418
airport_dist 4.0898 1 0.043
tourism 0.4325 1 0.511
healthcare 1.0549 1 0.304
popdens 0.0785 1 0.779
GLOBAL 5.3354 5 0.376
chisq df p
pers 0.5507 1 0.4580
age 3.8479 1 0.0498
male 5.6687 1 0.0173
conservative 2.3884 1 0.1222
academics 0.0979 1 0.7544
medinc 3.2705 1 0.0705
manufact 2.2318 1 0.1352
airport_dist 6.7326 1 0.0095
tourism 1.5199 1 0.2176
healthcare 1.5803 1 0.2087
popdens 0.0102 1 0.9194
onset_prev 3.0841 1 0.0791
slope_prev 4.5356 1 0.0332
GLOBAL 18.3991 13 0.1429
Assumptions US socdist onsets
list_iterater(us_list_results$us_cox_socdist_cpt, test = 'ph')
chisq df p
pers 40.1 1 2.4e-10
GLOBAL 40.1 1 2.4e-10
chisq df p
pers 40.37 1 2.1e-10
age 3.15 1 0.07575
male 8.88 1 0.00288
conservative 12.21 1 0.00047
GLOBAL 47.17 4 1.4e-09
chisq df p
pers 35.276 1 2.9e-09
academics 10.776 1 0.001
medinc 0.712 1 0.399
manufact 6.477 1 0.011
GLOBAL 37.669 4 1.3e-07
chisq df p
pers 45.88 1 1.3e-11
airport_dist 33.83 1 6.0e-09
tourism 25.00 1 5.7e-07
healthcare 1.81 1 0.18
popdens 47.05 1 6.9e-12
GLOBAL 101.93 5 < 2e-16
chisq df p
pers 34.785 1 3.7e-09
age 1.678 1 0.1953
male 6.300 1 0.0121
conservative 15.885 1 6.7e-05
academics 10.398 1 0.0013
medinc 0.941 1 0.3321
manufact 7.793 1 0.0052
airport_dist 29.969 1 4.4e-08
tourism 19.239 1 1.2e-05
healthcare 0.718 1 0.3968
popdens 34.662 1 3.9e-09
onset_prev 52.981 1 3.4e-13
slope_prev 59.230 1 1.4e-14
GLOBAL 121.454 13 < 2e-16
chisq df p
pers 12.4 1 0.00043
GLOBAL 12.4 1 0.00043
chisq df p
pers 12.03 1 0.00052
age 2.36 1 0.12437
male 6.84 1 0.00890
conservative 10.00 1 0.00157
GLOBAL 26.40 4 2.6e-05
chisq df p
pers 10.018 1 0.0016
academics 8.193 1 0.0042
medinc 0.203 1 0.6524
manufact 5.864 1 0.0155
GLOBAL 26.466 4 2.5e-05
chisq df p
pers 8.91 1 0.0028
airport_dist 35.75 1 2.2e-09
tourism 24.41 1 7.8e-07
healthcare 1.13 1 0.2872
popdens 63.76 1 1.4e-15
GLOBAL 112.08 5 < 2e-16
chisq df p
pers 8.922 1 0.00282
age 1.096 1 0.29514
male 5.384 1 0.02032
conservative 14.733 1 0.00012
academics 8.988 1 0.00272
medinc 0.515 1 0.47281
manufact 8.175 1 0.00425
airport_dist 32.147 1 1.4e-08
tourism 18.978 1 1.3e-05
healthcare 0.414 1 0.51984
popdens 36.422 1 1.6e-09
onset_prev 49.691 1 1.8e-12
slope_prev 59.684 1 1.1e-14
GLOBAL 120.284 13 < 2e-16
chisq df p
pers 0.561 1 0.45
GLOBAL 0.561 1 0.45
chisq df p
pers 0.427 1 0.51355
age 2.983 1 0.08417
male 7.640 1 0.00571
conservative 10.147 1 0.00145
GLOBAL 18.922 4 0.00081
chisq df p
pers 0.0957 1 0.7570
academics 9.7760 1 0.0018
medinc 0.5531 1 0.4571
manufact 6.6345 1 0.0100
GLOBAL 17.1795 4 0.0018
chisq df p
pers 0.402 1 0.53
airport_dist 32.089 1 1.5e-08
tourism 26.039 1 3.3e-07
healthcare 1.454 1 0.23
popdens 59.669 1 1.1e-14
GLOBAL 102.636 5 < 2e-16
chisq df p
pers 1.28e-03 1 0.97144
age 1.39e+00 1 0.23869
male 5.95e+00 1 0.01473
conservative 1.48e+01 1 0.00012
academics 1.03e+01 1 0.00133
medinc 8.99e-01 1 0.34313
manufact 8.21e+00 1 0.00417
airport_dist 2.96e+01 1 5.2e-08
tourism 1.96e+01 1 9.6e-06
healthcare 5.16e-01 1 0.47244
popdens 3.61e+01 1 1.9e-09
onset_prev 5.21e+01 1 5.2e-13
slope_prev 6.00e+01 1 9.3e-15
GLOBAL 1.18e+02 13 < 2e-16
chisq df p
pers 16.4 1 5.1e-05
GLOBAL 16.4 1 5.1e-05
chisq df p
pers 15.94 1 6.5e-05
age 2.83 1 0.0925
male 7.42 1 0.0065
conservative 10.10 1 0.0015
GLOBAL 28.85 4 8.4e-06
chisq df p
pers 13.905 1 0.00019
academics 9.124 1 0.00252
medinc 0.395 1 0.52976
manufact 6.149 1 0.01315
GLOBAL 32.382 4 1.6e-06
chisq df p
pers 11.5 1 0.0007
airport_dist 33.8 1 6.0e-09
tourism 24.3 1 8.3e-07
healthcare 1.1 1 0.2943
popdens 61.5 1 4.4e-15
GLOBAL 116.8 5 < 2e-16
chisq df p
pers 11.805 1 0.00059
age 1.281 1 0.25765
male 5.737 1 0.01661
conservative 14.695 1 0.00013
academics 9.849 1 0.00170
medinc 0.717 1 0.39707
manufact 8.389 1 0.00378
airport_dist 30.178 1 3.9e-08
tourism 18.997 1 1.3e-05
healthcare 0.393 1 0.53092
popdens 35.993 1 2.0e-09
onset_prev 51.526 1 7.1e-13
slope_prev 60.946 1 5.9e-15
GLOBAL 122.473 13 < 2e-16
chisq df p
pers 8.41 1 0.0037
GLOBAL 8.41 1 0.0037
chisq df p
pers 8.65 1 0.00328
age 1.99 1 0.15812
male 7.00 1 0.00814
conservative 10.05 1 0.00152
GLOBAL 19.68 4 0.00058
chisq df p
pers 9.659 1 0.00188
academics 7.704 1 0.00551
medinc 0.124 1 0.72482
manufact 7.032 1 0.00801
GLOBAL 20.333 4 0.00043
chisq df p
pers 7.99 1 0.0047
airport_dist 34.83 1 3.6e-09
tourism 24.38 1 7.9e-07
healthcare 0.57 1 0.4502
popdens 51.22 1 8.2e-13
GLOBAL 97.26 5 < 2e-16
chisq df p
pers 9.841 1 0.0017
age 0.987 1 0.3204
male 5.645 1 0.0175
conservative 15.669 1 7.5e-05
academics 9.134 1 0.0025
medinc 0.439 1 0.5078
manufact 9.371 1 0.0022
airport_dist 30.646 1 3.1e-08
tourism 19.652 1 9.3e-06
healthcare 0.180 1 0.6712
popdens 35.728 1 2.3e-09
onset_prev 50.955 1 9.4e-13
slope_prev 60.718 1 6.6e-15
GLOBAL 122.278 13 < 2e-16
Assumptions GER socdist adjustment levels
list_iterater(ger_list_results$ger_lm_socdist_mean, test = 'qq')

























list_iterater(ger_list_results$ger_lm_socdist_mean, test = 'bp')
studentized Breusch-Pagan test
data: .
BP = 6.6073, df = 1, p-value = 0.01016
studentized Breusch-Pagan test
data: .
BP = 11.402, df = 4, p-value = 0.02239
studentized Breusch-Pagan test
data: .
BP = 38.391, df = 4, p-value = 9.306e-08
studentized Breusch-Pagan test
data: .
BP = 0.91771, df = 5, p-value = 0.9689
studentized Breusch-Pagan test
data: .
BP = 36.57, df = 13, p-value = 0.0004836
studentized Breusch-Pagan test
data: .
BP = 2.8103, df = 1, p-value = 0.09366
studentized Breusch-Pagan test
data: .
BP = 16.321, df = 4, p-value = 0.002617
studentized Breusch-Pagan test
data: .
BP = 43.588, df = 4, p-value = 7.812e-09
studentized Breusch-Pagan test
data: .
BP = 3.9413, df = 5, p-value = 0.5579
studentized Breusch-Pagan test
data: .
BP = 40.982, df = 13, p-value = 9.594e-05
studentized Breusch-Pagan test
data: .
BP = 3.9066, df = 1, p-value = 0.0481
studentized Breusch-Pagan test
data: .
BP = 14.173, df = 4, p-value = 0.006764
studentized Breusch-Pagan test
data: .
BP = 37.668, df = 4, p-value = 1.312e-07
studentized Breusch-Pagan test
data: .
BP = 1.6777, df = 5, p-value = 0.8917
studentized Breusch-Pagan test
data: .
BP = 36.106, df = 13, p-value = 0.0005712
studentized Breusch-Pagan test
data: .
BP = 6.5055, df = 1, p-value = 0.01075
studentized Breusch-Pagan test
data: .
BP = 15.626, df = 4, p-value = 0.003565
studentized Breusch-Pagan test
data: .
BP = 40.792, df = 4, p-value = 2.968e-08
studentized Breusch-Pagan test
data: .
BP = 1.6497, df = 5, p-value = 0.8952
studentized Breusch-Pagan test
data: .
BP = 39.059, df = 13, p-value = 0.0001956
studentized Breusch-Pagan test
data: .
BP = 2.3883, df = 1, p-value = 0.1222
studentized Breusch-Pagan test
data: .
BP = 14.994, df = 4, p-value = 0.004713
studentized Breusch-Pagan test
data: .
BP = 38.685, df = 4, p-value = 8.091e-08
studentized Breusch-Pagan test
data: .
BP = 2.6498, df = 5, p-value = 0.7538
studentized Breusch-Pagan test
data: .
BP = 38.358, df = 13, p-value = 0.000253
Assumptions US socdist adjustment levels
list_iterater(us_list_results$us_lm_socdist_mean, test = 'qq')

























list_iterater(us_list_results$us_lm_socdist_mean, test = 'bp')
studentized Breusch-Pagan test
data: .
BP = 11.773, df = 1, p-value = 0.0006008
studentized Breusch-Pagan test
data: .
BP = 11.552, df = 4, p-value = 0.02101
studentized Breusch-Pagan test
data: .
BP = 9.536, df = 4, p-value = 0.04901
studentized Breusch-Pagan test
data: .
BP = 37.618, df = 5, p-value = 4.501e-07
studentized Breusch-Pagan test
data: .
BP = 33.974, df = 13, p-value = 0.001215
studentized Breusch-Pagan test
data: .
BP = 6.7876, df = 1, p-value = 0.00918
studentized Breusch-Pagan test
data: .
BP = 10.068, df = 4, p-value = 0.03929
studentized Breusch-Pagan test
data: .
BP = 11.336, df = 4, p-value = 0.02303
studentized Breusch-Pagan test
data: .
BP = 49.661, df = 5, p-value = 1.626e-09
studentized Breusch-Pagan test
data: .
BP = 34.866, df = 13, p-value = 0.0008878
studentized Breusch-Pagan test
data: .
BP = 0.18699, df = 1, p-value = 0.6654
studentized Breusch-Pagan test
data: .
BP = 14.373, df = 4, p-value = 0.006195
studentized Breusch-Pagan test
data: .
BP = 10.332, df = 4, p-value = 0.03519
studentized Breusch-Pagan test
data: .
BP = 50.528, df = 5, p-value = 1.081e-09
studentized Breusch-Pagan test
data: .
BP = 34.641, df = 13, p-value = 0.0009611
studentized Breusch-Pagan test
data: .
BP = 14.662, df = 1, p-value = 0.0001286
studentized Breusch-Pagan test
data: .
BP = 12.049, df = 4, p-value = 0.01699
studentized Breusch-Pagan test
data: .
BP = 15.089, df = 4, p-value = 0.004519
studentized Breusch-Pagan test
data: .
BP = 54.746, df = 5, p-value = 1.472e-10
studentized Breusch-Pagan test
data: .
BP = 36.335, df = 13, p-value = 0.0005262
studentized Breusch-Pagan test
data: .
BP = 9.3516, df = 1, p-value = 0.002228
studentized Breusch-Pagan test
data: .
BP = 20.051, df = 4, p-value = 0.000488
studentized Breusch-Pagan test
data: .
BP = 11.455, df = 4, p-value = 0.0219
studentized Breusch-Pagan test
data: .
BP = 51.643, df = 5, p-value = 6.386e-10
studentized Breusch-Pagan test
data: .
BP = 36.846, df = 13, p-value = 0.000438
LS0tCnRpdGxlOiAiQXNzdW1wdGlvbiBDaGVja3MiCm91dHB1dDogaHRtbF9ub3RlYm9vawotLS0KCgojIEltcG9ydCBtb2RlbCBvYmplY3RzIApgYGB7cn0KCmxvYWQoJy9Vc2Vycy9ocDI1MDAvR29vZ2xlIERyaXZlL1NUVURZL0NvbHVtYmlhL1Jlc2VhcmNoL0Nvcm9uYS9EYXRhL0dFUi9nZXJfbGlzdF9yZXN1bHRzX2ZpeGVkX3dpbmRvdy5SRGF0YScpCmxvYWQoJy9Vc2Vycy9ocDI1MDAvR29vZ2xlIERyaXZlL1NUVURZL0NvbHVtYmlhL1Jlc2VhcmNoL0Nvcm9uYS9EYXRhL1VTL3VzX2xpc3RfcmVzdWx0c19maXhlZF93aW5kb3cuUkRhdGEnKQoKbGlicmFyeShsbXRlc3QpCmxpYnJhcnkoY2FyKQpsaWJyYXJ5KHN1cnZpdmFsKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKCgpgYGAKCiMgRGVmaW5lIEZ1bmN0aW9ucyAKYGBge3J9CgpsaXN0X2l0ZXJhdGVyIDwtIGZ1bmN0aW9uKG1vZGVscywgdGVzdCkgewogIGZvcihpIGluIG1vZGVscyl7CiAgICBmb3IoaiBpbiBpKXsKICAgICAgCiAgICAgIGlmKHRlc3QgPT0gJ3FxJyl7aiAlPiUgcGxvdCgyKX0KICAgICAgaWYodGVzdCA9PSAna3MnKXtqICU+JSByZXNpZCgpICU+JSBrcy50ZXN0KHk9cG5vcm0pICU+JSBwcmludCgpfQogICAgICBpZih0ZXN0ID09ICdicCcpe2ogJT4lIGJwdGVzdCgpICU+JSBwcmludCgpfQogICAgICBpZih0ZXN0ID09ICdwaCcpe2ogJT4lIGNveC56cGgoKSAlPiUgcHJpbnQoKX0KICAgIH0KICB9Cn0KYGBgCgoKIyBBc3N1bXB0aW9ucyBHRVIgQ09WSUQtMTkgb25zZXRzIChwcm9wb3J0aW9uYWwgaGF6YXJkcykKYGBge3J9Cmxpc3RfaXRlcmF0ZXIoZ2VyX2xpc3RfcmVzdWx0cyRnZXJfY294X3ByZXZfb25zZXQsIHRlc3QgPSAncGgnKQpgYGAKCiMgQXNzdW1wdGlvbnMgVVMgQ09WSUQtMTkgb25zZXRzIChwcm9wb3J0aW9uYWwgaGF6YXJkcykKYGBge3J9Cmxpc3RfaXRlcmF0ZXIodXNfbGlzdF9yZXN1bHRzJHVzX2NveF9wcmV2X29uc2V0LCB0ZXN0ID0gJ3BoJykKYGBgCgojIEFzc3VtcHRpb25zIEdFUiBDT1ZJRC0xOSBncm93dGggcmF0ZXMgKG5vcm1hbGl0eSBvZiByZXNpZHVhbHMpCmBgYHtyfQpsaXN0X2l0ZXJhdGVyKGdlcl9saXN0X3Jlc3VsdHMkZ2VyX2xtX3ByZXZfc2xvcGUsIHRlc3QgPSAncXEnKQpgYGAKCmBgYHtyfQpsaXN0X2l0ZXJhdGVyKGdlcl9saXN0X3Jlc3VsdHMkZ2VyX2xtX3ByZXZfc2xvcGUsIHRlc3QgPSAnYnAnKQpgYGAKCiMgQXNzdW1wdGlvbnMgVVMgQ09WSUQtMTkgZ3Jvd3RoIHJhdGVzIChub3JtYWxpdHkgb2YgcmVzaWR1YWxzKQpgYGB7cn0KbGlzdF9pdGVyYXRlcih1c19saXN0X3Jlc3VsdHMkdXNfbG1fcHJldl9zbG9wZSwgdGVzdCA9ICdxcScpCmBgYAoKYGBge3J9Cmxpc3RfaXRlcmF0ZXIodXNfbGlzdF9yZXN1bHRzJHVzX2xtX3ByZXZfc2xvcGUsIHRlc3QgPSAnYnAnKQpgYGAKCiMgQXNzdW1wdGlvbnMgR0VSIHNvY2Rpc3Qgb25zZXRzCmBgYHtyfQpsaXN0X2l0ZXJhdGVyKGdlcl9saXN0X3Jlc3VsdHMkZ2VyX2NveF9zb2NkaXN0X2NwdCwgdGVzdCA9ICdwaCcpCmBgYAoKIyBBc3N1bXB0aW9ucyBVUyBzb2NkaXN0IG9uc2V0cwpgYGB7cn0KbGlzdF9pdGVyYXRlcih1c19saXN0X3Jlc3VsdHMkdXNfY294X3NvY2Rpc3RfY3B0LCB0ZXN0ID0gJ3BoJykKYGBgCgojIEFzc3VtcHRpb25zIEdFUiBzb2NkaXN0IGFkanVzdG1lbnQgbGV2ZWxzCmBgYHtyfQpsaXN0X2l0ZXJhdGVyKGdlcl9saXN0X3Jlc3VsdHMkZ2VyX2xtX3NvY2Rpc3RfbWVhbiwgdGVzdCA9ICdxcScpCmBgYAoKYGBge3J9Cmxpc3RfaXRlcmF0ZXIoZ2VyX2xpc3RfcmVzdWx0cyRnZXJfbG1fc29jZGlzdF9tZWFuLCB0ZXN0ID0gJ2JwJykKYGBgCgojIEFzc3VtcHRpb25zIFVTIHNvY2Rpc3QgYWRqdXN0bWVudCBsZXZlbHMKYGBge3J9Cmxpc3RfaXRlcmF0ZXIodXNfbGlzdF9yZXN1bHRzJHVzX2xtX3NvY2Rpc3RfbWVhbiwgdGVzdCA9ICdxcScpCmBgYAoKYGBge3J9Cmxpc3RfaXRlcmF0ZXIodXNfbGlzdF9yZXN1bHRzJHVzX2xtX3NvY2Rpc3RfbWVhbiwgdGVzdCA9ICdicCcpCmBgYAo=